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We report the first observation of Xci""*, a Cabibbo- and color-suppressed decay in a data 

sample of 386 x 10® BB events collected at the T(4S') resonance with the Belle detector at the KEKB 
asymmetric-energy e"'"e~ collider. We observe 55 ± 10 signal events with a statistical significance of 
6.3(T including systematic uncertainties. The measured branching fraction and charge-asymmetry 
is XciTT"^) = (2.2 ± 0.4 ± 0.3) x lO"'^ and A-^ = 0.07 ± 0.18 ± 0.02, respectively. We also 

determine the ratio ^ Xci7r=^)/S(B± ^ XaK"^) = (4.3 ± 0.8 ± 0.3)%. 

PACS numbers: 13.25.Hw, 14.40Gx, 14.40.Nd 



Decays of iJ-mesons to two-body final states including 
charmonium are expected to occur predominantly via the 
color-suppressed spectator diagram as shown in Fig. 1. 
The branching fraction for the B~ XciK~ decay 
mode is well measured by Belle and BaBar [^|^. To pro- 
duce this final state, the vector current {W~) couples to 
a cs pair; the s-quark and spectator anti-quark hadronize 
into a kaon. If this theoretical description is correct, 
a corresponding Cabibbo-suppressed decay mode should 
exist, where the vector current couples to a cd pair. The 
d-quark and spectator anti-quark hadronize as a pion, 
which leads to a B~ Xci'^~ decay. If the leading- 
order tree level diagram is the dominant contribution, 
the factorization picture implies that the branching frac- 
tion of B~ Xcii"" decay mode should be ~ 5% of that 
of the Cabibbo-allowed B~ XciK~ decay mode 0. 
The Standard Model predicts that for b ccs decays, 
the tree and penguin contributions have a small relative 
weak phase. Therefore, negligible direct CP- violation is 
expected in B^ — > XciK^ decay. In 6 ^ ccd transitions, 
however, tree and penguin contributions have different 
phases and direct CP-violation may be as large as a few 
percent 0, |^ . 

In this paper, we report the first observation of B~ — s- 
XciTT^ decay. A measurement of the ratio of branch- 
ing fractions B{B~ Xci'^~)/^{B^ xaK^) and 
a search for direct CP-violation in B^ Xci""* de- 
cays is also presented. We use a data sample con- 
taining (386 ± 5) X 10^ BB events collected at the 
T(4S') resonance with the Belle detector at the KEKB 
asymmetric-energy e+e^ collider Q. 

The Belle detector is a large solid- angle magnetic spec- 
trometer located at the KEKB e+e^ storage rings, which 
collide 8.0 GeV electrons with 3.5 GcV positrons pro- 
ducing a center-of-mass (CM) energy of 10.58 GeV, the 
mass of the T{AS) resonance. Closest to the interaction 
point (IP) is a silicon vertex detector (SVD), surrounded 
by a 50-layer central drift chamber (CDC), an array 
of aerogel Cherenkov counters (ACC), a barrel- like ar- 




FIG. 1: Leading-order tree level diagram for the decays under 
study. 



rangement of time-of-flight (TOF) scintillation counters, 
and an electromagnetic calorimeter (ECL) comprised of 
CsI(T^) crystals. These subdetectors are located inside a 
superconducting solenoid coil that provides a 1.5 T mag- 
netic field. An iron flux-return yoke located outside the 
coil is instrumented to detect mesons and to identify 
muons. The detector is described in detail elsewhere 
The data set consists of two subsets: the first 152 x 10^ 
P-meson pairs were collected with a 2.0 cm radius beam- 
pipe and a 3-layer SVD, and the remaining 234 x 10^ 
P-meson pairs with a 1.5 cm radius beam-pipe, a 4-layer 
SVD and a small-cell inner drift chamber j9j. 

Events with P-meson candidates are first selected by 
applying general hadronic event selection criteria. These 
include a requirement on charged tracks (at least three 
of them should originate from an event vertex consistent 
with the IP), a requirement on the reconstructed CM en- 
ergy [E^^^ > 0.2v^, where ^/s is the total CM energy), a 
requirement on the longitudinal (z-direction) component 
of the reconstructed CM momentum with respect to the 
beam direction (|p^^^| < Q.5y/s/c), and a requirement on 
the total ECL energy (O.l^i < Pg^i < 0.8Vi) with at 
least two energy clusters. To suppress continuum back- 
ground, we reject events where the ratio of the second to 
zeroth Fox- Wolfram moments is greater than 0.5. To 
remove charged particle tracks that are poorly measured 
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FIG. 2: The AM {Mi,+i-^ - distribution for the Xci 

candidates. The arrows indicate the selected mass region. 
The enhancement just above the Xci mass region is due to 
the Xc2. 



or do not come from the interaction region, we require 
their origin to be within 0.5 cm of the IP in the radial di- 
rection, and 5 cm along the beam direction (z-direction). 

We reconstruct the Xci state via the decay mode 
Xci li>- We begin by reconstructing J — > 

candidates, where ^ is a muon or electron. For muon 
tracks, identification is based on track penetration depth 
and the hit pattern in the KLM system. Electron tracks 
are identified by a combination of dE/dx from the CDC, 
E/p (E is the energy deposited in the ECL and p is 
the momentum measured by the SVD and the CDC), 
and shower shape in the ECL. In order to recover di- 
electron events in which one or both electrons radiate 
a photon, the four-momenta of all photons within 0.05 
radians of the e+ or e~ directions are included in the 
invariant mass calculation. The invariant mass win- 
dow used to select J/4> candidates in the /.t^/i~(e+e^) 
channel is -0.06 (-0.15) GeV/c^ < M^+f- - rrij/^ < 
0.036 GeV/c^, where mj/^ denotes the nominal J/ip 
mass llj; these intervals are asymmetric in order to 
include part of the radiative tails. Vertex- and mass- 
constrained kinematic fits are performed for selected J/^p 
candidates to improve the momentum resolution. 

Photons are identified as ECL energy clusters that are 
not associated with a charged track and have a minimum 
energy of 0.060 GeV. We reject the photon candidate if 
the ratio of the energy in the array of the central 3x3 
ECL cells to that in the array of 5 x 5 cells is less than 
0.87. 

To reconstruct the Xci state, we combine a J/^ candi- 
date with momentum below 2.0 GeV/c in the CM frame 
with a selected photon. To suppress photons arising from 
7r° 77, we veto photons that, when combined with an- 



other photon in the event, satisfy 0.110 GeV/c^ < M^^ < 
0.150 GeV/c^. The Xci candidates are selected by requir- 
ing the mass difference (AM = Mg+i-^ — Mg+i- ) to he 
between 0.370 GeV/c^ and 0.438 GeV/c^. The AM dis- 
tribution is shown in Fig. 2. A mass-constrained fit is 
applied to all selected Xci candidates in order to improve 
the momentum resolution. 

Charged pions and kaons are identified using energy 
loss measurements in the CDC, Cherenkov light yields in 
the ACC, and TOF information. The information from 
these detectors is combined to form n-K likelihood ratio, 
TZ{tt/K) ~ Lt^I{Lt^ -I- Ck), where £^(£if) is the likeli- 
hood that a pion (kaon) would produce the observed de- 
tector response. Charged tracks with TZ{'k/K) > 0.9 are 
selected as charged pions, and tracks with TZ{tt/K) < 0.4 
are selected as charged kaons. The efficiency for pion 
(kaon) identification is 75.1% (86.1%) and the probabil- 
ity of kaon (pion) misidentification is 4.6% (10.5%) with 
the above criteria. We determine the selection criteria by 
optimizing the figure of merit, S/^/{S + B), where S{B) 
is the number of signal (background) events in the signal 
region, with an assumed branching fraction that is 5% of 
that for B- -> XciK' [ij. 

We reconstruct i3-mesons by combining a Xci can- 
didate with a charged pion or kaon. The energy dif- 
ference, AE' = — i?beam ^^"^ ^^'^ bcam-coustrained 
mass Mbc = \/ -E-bcam ~ P*B ' ^'^^ v^sed to separate sig- 
nal from background, where E^^^^ is the run dependent 
beam energy, and E'^ and p*^ are the reconstructed en- 
ergy and momentum, respectively of the B-meson can- 
didates in the CM frame. We accept candidates in 
the region 5.27 GeV/c^ < Mbc < 5.29 GeV/c^ and 
|A£;| < 0.2 (0.15) GeV for the B' Xci'^~{K-) mode. 
When an event contains more than one B-meson candi- 
date passing the above requirements (this occurs in ~ 
2.5% of the candidate events), the candidate with Afbc 
closest to the nominal B~ mass is selected. 

We extract the signal yields by performing a binned 
maximum likelihood fit to the Ai? distribution of the se- 
lected candidates. For the B^ XciK^ mode, we fit 
with a sum of two Gaussians for signal and a second- 
order polynomial for background. In the fit for the 
B^ — > XciTT" mode, a background component exists due 
to misidentified Xci^^~. This background has a peak 
at AE —0.07 GeV and is modeled by a sum of two 
bifurcated Gaussians. A third-order polynomial is used 
for the sum of all other backgrounds. We study back- 
grounds using a large sample of inclusive charmonium 
Monte Carlo (MC) events TZ|. Except for the misidenti- 
fied B~ — > XciK~ background, no structure is observed 
in the AE distribution. However, the Mbc distribution 
has a peaking component. Therefore, we use the AE dis- 
tributions for the signal extraction. The scatter plot of 
AE versus Mbc for B^ Xcitt^ candidates is shown 
in Fig. 3, where the Afbc requirement is loosened to 
5.2 GeV/c2. 

All parameters of the fitting functions are fioated in 
the fit of the B^ XciK~ mode. For the B~ Xci'^~ 
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FIG. 3: The scatter plot of versus Mbc for Xcitt''' 
candidates, where the two vertical lines indicate the Afbc re- 
gion used for signal extraction. 



FIG. 4: The AE distribution for the XaK"^ decay 

mode. The solid and dashed curves show the total fit and the 
polynomial background component of the fit, respectively. 



mode, the signal shape is fixed to that obtained from the 
XciK^ mode. The shape of misidentified XciK~' 
background is initially determined from a MC sample, 
with a correction applied to account for the small differ- 
ence between data and MC in the XciK^ sam- 
ple. We obtain 1597 ± 48 and 55 ± 10 signal events for 
B~ — !■ XciK~ and B~ Xcitt" modes, respectively. 
The AE distributions are shown in Figs. 4 and 5, to- 
gether with the fit results. The number of of misidenti- 
fied B^ XciK^ events obtained from the fit to Fig. 5 
is 61 ± 14. This is consistent with the expectation from 
the observed B~ XciK~ signal yield (Fig. 4) given the 
probability of misidentifying a kaon as a pion. The signif- 
icance of the B~ — > Xci'"'" signal is 6.3ct, where the signif- 
icance is defined as a/— 21n(£o/'Cmax) and £max ('Co) de- 
notes the likelihood value at the maximum (with the sig- 
nal yield fixed at zero) . We include the effect of system- 
atic error in this calculation by subtracting a quadratic 
sum of the variations of the significance in smaller direc- 
tion when each fixed parameter in the fit is changed by 
±lcr. 

The branching fraction for the B~ Xcit^~ decay 
mode is calculated by dividing the observed signal yield 
by the reconstruction efhciency, the number of BB events 
in the data sample, and the daughter branching fractions. 
We determine the reconstruction efficiency (17.3%) from 
signal MC events, where the correction for difference be- 
tween data and MC has been applied for the pion identifi- 
cation requirement (0.90±0.01) and the AM requirement 
(0.97±0.03). The correction factor for the AM require- 
ment is determined from the B~ — s- XciK~ sample and 
is estimated by taking the ratio of yields from data and 
MC for tight (0.370 GeV/c^ < AM < 0.438 GeV/c^) 
and loose (0.3 GeV/c^ < AM < 0.5 GeV/c^) AM win- 
dows. We use the daughter branching fractions published 
in Particle Data Book 2004 Equal production of 

neutral and charged i?-meson pairs in T(4S) decay is as- 
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FIG. 5: The AE distribution for the Xc-i^^^ decay 

mode. The signal peak is seen around zero. The peak at 
—0.07 GeV is from B^ XciK^ decay. The solid and dashed 
curves show the total fit and the polynomial background com- 
ponent of the fit, respectively. 



sumed. The resulting branching fraction is 

B{B^ ^ XciTT^) = (2.2 ± 0.4 ± 0.3) x 10"^ (1) 

where the first error is statistical and the second is 
systematic. We obtain the branching fraction for the 
B~ XciK~ decay mode by similar procedures. The 
result, (51.4 ± 1.5) x 10~^, (error is statistical only) is 
consistent with the previous measurements 0, 1^ ■ The 
ratio of branching fractions is 



BjB- ^ XciTT-) 
B{B- ^ XciK-) 



(4.3 ±0.8 ±0.3)%, 



(2) 



which is consistent with expectations from the factoriza- 
tion model Q. 

The systematic uncertainties are summarized in Ta- 
ble ^ Since the shape used for the signal is fixed in the 
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TABLE I: Summary of systematic errors on branching frac- 
tion. 



Source 



Uncertainty (%) 



Uncertainty in yield 
Tracking error 
Lepton Identification 
PID (pion) 
7 detection 
AM requirement 
MC statistics 
Nbb 

Daughter branching fractions 



5.9 
3.0 
4.0 
1.0 
2.0 
3.0 
0.9 
1.2 
10.6 



yield extraction, 0.007; possible difference between B~ 
and 5+ signal shape parameters, 0.002 (0.001); possi- 
ble charge asymmetry in pion (kaon) identification effi- 
ciency 0.014 (0.011); and possible detector bias 0.016, 
which is estimated from the charge asymmetry of the 
J/^pK^ decay sample without a PID requirement. 



TABLE II: Summary of charge asymmetries. First (second) 
error is statistical (systematic). 



Mode 



Total 



13.7 



Yield(-) Yield(+) 



A 



29 ±7 
792 ± 31 



25 ±7 
807 ± 31 



0.07 ±0.18 ±0.02 
-0.01 ±0.03 ±0.02 



fit to Fig. 5, we repeat the fit, varying each fixed shape 
parameter by the ilcr uncertainty in our determination 
of it from external samples (Fig. 4 and MC). The system- 
atic uncertainty on the signal yield is then calculated by 
taking the quadratic sum of the deviations in the signal 
yield from the nominal value. We checked for possible 
bias in the fitting using a MC sample; no significant bias 
was found. The systematic uncertainty assigned to the 
yield is 5.9%. The uncertainty on the tracking efficiency 
is estimated to be 1.0% per track, while that due to lep- 
ton identification is 2.0% per lepton, and 1.0% (1.3%) 
per pion (kaon) identification (PID). We assign an un- 
certainty of 2.0% for the 7 detection efficiency. The sys- 
tematic uncertainty due to Xci iJ 14' and J/ij} — > 
branching fractions is 10.6%. The total systematic error 
is the sum of all the above uncertainties in quadrature. 

Many of the systematic errors cancel for the ratio of 
branching fractions; contributions come from only the 
uncertainty in the Xcitt^ yield, PID (1.0% for 

B- XciK^), and MC statistics (1.0% for B' 
XciK-). 

The CP- violating charge asymmetry Ai is defined as 



N7 



i = 77, K. 



(3) 



Here, N~ and Nl^ are the signal yields for negative and 
positive i?-meson decays and are measured separately by 
using the method described above. For the B^ Xci^r" 
mode, the polynomial background shape is fixed to that 
obtained for the branching fraction measurement. 

The measured charge asymmetries for the B^ — s- 
Xci7r^(ivr*) decay modes are listed in Table Hll No signif- 
icant asymmetries are seen in either decay modes. The 
systematic errors on A^r (Ak) include: uncertainty in 



In summary, we report the first observation of B^ — > 
XciTT" decay with 386 x 10^ BB events. The observed sig- 
nal yield is 55 ± 10 with a significance of 6.3ct including 
systematic uncertainty. The measured branching fraction 



is S(B^ 



XclTT-' 



(2.2 ± 0.4 ± 0.3) X 10-5. The ratio 



B{B- xci^-) IB{B- ^ XciK-) = (4.3 ±0.8 ±0.3)%, 
which is consistent with the Standard Model prediction. 
While the accuracy of Ak is improved from the previ- 
ous measurement no significant CP-violating charge 
asymmetries are observed in either B^ — > Xci""^ or 
B^ XciK"^ decay modes. 
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